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ABSTRACT 



A gas -cooled, cylindrical, nuclear reactor is used as the 
basis in the theoretical derivation of the coolant's enthalpy increase 
and pressure loss across the reactor. Turbulent flow is assumed in 
the coolant passages of the reactor, and the Reynolds analogy is 
used as the basis for correlating the heat transfer coefficient and 
the friction coefficient. 

The general equations are derived and two examples, a 
nuclear -hydrogen rocket and a ramjet, are given to demonstrate 
applications of the general results. 



SUMMARY 



The development of nuclear reactors has led to a new class of heat 
transfer problems. Heat is produced within a solid fuel at a rate that 
is almost independent of the coolant flow rate, at least for gas -cooled 
reactors. The power density varies with position in general, although 
it can be made nearly uniform with the introduction of appropriate re- 
flectors and with variation of fuel concentration. However, these 
measures are not desirable, particularly if there are weight restric- 
tions for the reflector and cost restrictions on fuel elements. 

Aside from considerations of weight and cost, the most important 
characteristic of a reactor is the ratio of pumping work to the increase 
of thermodynamic availability of the coolant passing through the re- 
actor. For any single coolant passage through the reactor, this ratio 
can be determined readily from the Reynolds analogy, for any speci- 
fied power distribution along the length of the channel, and a given 
coolant flow rate through the channel or a given pressure drop over the 
reactor. The limitation on the solid fuel temperature will fix the 
minimum pressure drop or flow rate. If the power density is the same 
for all channels, the results are simply summed over the reactor 
cross-section. 

When the power density varies from channel to channel, there is 
no simple way of calculating pressure drop or temperature rise 
through the reactor. The pressure drop for all channels must be the 
same, but the flow rate and temperature rise will vary from channel to 
channel. Apparently the calculations of temperature and flow rate dis- 
tributions over the reactor have been done by numerical computation 
in the past. An iterative procedure is developed in this thesis for 



determination of the overall reactor performance in analytical form. 

The sole restriction on the procedure ie that the Mach number of the 
flow entering any coolant channel be appreciably less than unity. This 
restriction is not a serious one because the pressure drop resulting 
from higher entering Mach numbers would be generally unacceptable. 

The principal objective of this thesis is the development of the 
method of calculating reactor performance for any specified power den- 
sity distribution. In final form the method gives the pressure drop, 
the temperature distribution and mass flow rate distribution in the 
coolant leaving the reactor, as well as maximum wall temperatures, 
where porosity, coolant channel geometry, and power density distribu- 
tion are prescribed. Tho temperature of the coolant leaving the re- 
actor can be made uniform over the reactor cross-section by introducing 
appropriate restrictions in the Individual coolant channels. This pro- 
cedure increases the average temperature rise through the reactor for 
a given maximum wall temperature, but does so at the cost of in- 
creased pressure drop. Whether restriction of the coolant flow channels 
in this way is advantageous or not will depend upon the particular ap- 
plication. The method of analysis developed here can be used to deter- 
mine the optimum restriction of the flow channels for any specific 
example. 

To illustrate the application of the method for calculating reactor 
performance, the nuclear rocket and nuclear ramjet were chosen as 
examples. The specific impulse for the rocket is found to be higher 
when the coolant passages are restricted to give a uniform coolant out- 
let temperature than for no restriction. The thrust coefficient of the 
ramjet ie not improved appreciably by restricting the flow rate. 
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SYMBOLS 



A 

a. 



c 

P 



m 



k 

/ 

rh 



p 



p 

q 

R 

T 



u 



cC 

P 

Pf 

7 

Z 4 

i) 



Coolant passage cross -sectional area 
Radius of reactor 
Friction coefficient 

Heat transfer coefficient (Stanton number) 

Specific heat at constant pressure 

Axial power distribution function 

Radial power distribution function 

Conductivity 

Coolant passage length 

Total mass flow rate through reactor 

Mass flow rate per unit area 

Coolant passage wetted perimeter 

Gas pressure (static) 

Power density in reactor; heat produced per unit volume 
Maximum power density in reactor 
Gas constant 

Gas temperature (static, averaged over coolant passage 
cross-section) 

Gas velocity (averaged over coolant passage cross- 
section) 

Porosity; void volume /total volume 

p ✓ 

C h ~£- , heat transfer parameter 
r 

f F/ , friction parameter 

T 7T 

r / cl , dimensionless radius from centerline of reactor 
Viscosity 

Kinematic viscosity 

x/Z , dimensionless distance from coolant passage 
entrance 



$ 



SYMBOLS 



p Gas density (averaged over coolant passage cross- 

section) 

^ Prandtl number 



Subscripts: 

0 

1 

w 

T 



Entrance to coolant passage 
Exit from coolant passage 
Wall of coolant passage 
Total 



I. INTRODUCTION 



The power distribution in a nuclear reactor is proportional to 
the neutron density which, in a bare cylindrical reactor, is highest in 
the center and decreases both axially and radially av/ay from the cen- 
ter. When the reactor is used as a source of heat, the variation in 
power density distribution results in an uneven heat supply to the 
coolant. 

The neutron density can be made more uniform by enclosing the 
reactor with a reflector which causes some of the leaking neutrons to 
diffuse back into the core, or by varying the fuel loading so that the 
fuel concentration is higher towards the edges than in the center. 
However, both of these methods have disadvantages. The reflector 
material is usually quite bulky and relatively heavy. As a result, in 
several applications of nuclear reactors as a source of power, spe- 
cifically, for airborne propulsive applications, it is necessary to 
restrict the amount of reflector material that might normally be used. 
Fuel concentration variation is undesirable since reactors are usually 
constructed by assembling many small lamina, or rods, to obtain the 
final configuration. For manufacturing ease and interchangeability, 
it is more desirable to have uniformity of the fuel elements. 

The uneven power density distribution of the nuclear reactor 
results In the coolant having temperature gradients in the exhaust if 
the mass flow in each of the coolant channels is the same. The cool- 
ant in the channel through the center of the core will be heated the 
most, whereas the flow in the outer channels will be much cooler due 
to the lower heat generation at the reactor extremities. Flow of this 
type is usually undesirable, since the surface temperature of the chan- 
nel wall 3 in the outer channels is much lower than the maximum 
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allowable, and thus the overall efficiency of the power plant is reduced. 
In addition, the temperature gradients of the exhaust of this type flow 
lead to mixing losses in the nozzle. 

Theoretically, these exhaust gas temperature gradients can be 
eliminated by restricting the mass flow in the outer channels so that 
the mass flow in every channel is proportional to the heat release in 
that channel. The result is that the outer channels, where the heat 
release is less, are cooled to a lesser degree so that each channel 
operates at its maximum allowable wall surface temperature, and the 
temperature of the exhaust gases from the reactor is uniform. This 
mass flow restriction can be accomplished by varying the flow channel 
cross-sectional areas continuously across the reactor core so that for 
the same overall pressure drop across the channel the mass flow is 
smellier. The same effect can be obtained by installing orifices in the 
channels to induce a pressure drop so that again the mass flow is pro- 
portional to the heat release of the particular channel. Of these two 
methods, the former is preferable since it utilizes all the available 
pressure drop to improve the heat transfer coefficient. However, it 
would not usually be practical to assemble a reactor in this manner, 
and so the orificing technique is considered to be more practical. 
Elimination of the exit temperature gradients and increasing the wall 
surface temperatures of all channels to the maximum allowable by 
either of these methods increases the overall pressure drop across 
the reactor and reduces the overall mass flow, in addition to increas- 
ing the average temperature rise per unit mass flow. Because of the 
increased pressure losses and the reduction in mass flow, the over- 
all performance of the reactor must be evaluated to determine the 



performance gain, if any, by orificing tc ascertain if the increased 
design and assembly complications are warranted. 

W. B. Hall (1) has worked out an example of the restricted flow 
reactor configuration. The basic reactor assumptions made by Hall 
and those of this thesis are similar. However, the analysis here is 
more general, including allowance for higher pressure drop and 
higher velocities in the coolant channels, and a procedure for finding 
overall performance when the coolant outlet temperature varies 
radially. These improvements in the analysis are needed for appli- 
cation to reactors for propulsion. 

In order to keep the general discussion of heat transfer and 

pressure drop in gas -cooled reactors for propulsive applications with- 

> 

in reasonable limits, simplifying assumptions are required. For 
this report, the reactor was taken a3 cylindrical, with radial sym- 
metry in fuel loading and coolant passage distribution. The coolant 
passages were assumed to be of uniform 3ize, and the cro83-sectional 
area was assumed constant over the length of the reactor. The power 
density of the reactor was assumed to vary axially and radially ac- 
cording to simple algebraic laws. Since the analysis was made pri- 
marily to Investigate the result of using a reactor as a type of heat 
exchanger, or source of heat, reactor physics as such was not con- 
sidered, and no account was taken of the power distribution shift 
during operation, etc. The porosity of the reactor was not assumed 
constant in the derivation of the general equations, but was assumed 
constant in specific examples worked out to demonstrate the use of 



the results. 



The evaluation of the overall performance of a propulsion 
system involves determining the overall enthalpy rise of the working 
fluid and the losses of the system. For a. reactor of given design, 
the controlling factor in operating power is normally some maximum 
allowable temperature in the system. This limit may be set by the 
phase stability of the fuel elements, or the moderator, by the allow- 
able thermal stresses in the fuel or some other part of the system, by 
the influence of temperature on corrosion, or some other particular 
thermal effect. For a given design and this limiting temperature, 
the temperature distribution in the reactor core can be determined, 
and from this, the maximum allowable coolant wall temperature 
subject to the limiting temperature of the particular configuration. 

Since this report was intended to be quite general, the controlling 
factor for reactor operation was selected as the maximum allowable 
coolant wall temperature which can be determined when given a par- 
ticular configuration with type of fuel, concentration, etc. Thus the 
problem to be solved involves calculating the overall enthalpy rise 
and preesure losses for a given maximum wall temperature and given 
initial conditions. 

The high power requirements of today's rockets and ramjets 
and the expanding technology in the nuclear reactor field make the pos- 
sibility of using the gas -cooled reactor as a source of propulsive 
power for these vehicles a distinct reality. Applications for rocket 
engines hold promise, since the effective exhaust velocity is propor- 
tional to the inverse of the square root of the molecular weight, which 
suggests the possibility of using low molecular weight gases ouch as 
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hydrogen for the working fluid. For ramjet applications, the nuclear 
reactor offers the possibility of relatively unlimited range or endur- 
ance at low as well ao high altitudes. 

In this thesis, the general equations are derived; and two exam- 
ples, a rocket using hydrogen ax the working fluid and a ramjet, are 
given. The results of these examples show a distinct advantage in ori- 
ficing to increase the specific impulse of a rocket, but the advantages 
for a ramjet are so slight that it is doubtful if the additional complexity 
in reactor construction arising from orificing is warranted. 
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H. ANALYSIS AND TENZE.AL SOLUTION 



2 . 1 Reyn olds An alog y 

In moat nuclear reactors the flow of the coolant through the 
cooling passages is turbulent. For this reason the motion of tho fluid 
cannot be defined exactly, 30 exact analytical solution of flow and heat 
transfer problems is impossible. However’, using the equations of 
motion, where the laminar transport properties are replaced by ef- 
fective value 8 which account for the effects of turbulence and which 
must be determined experimentally, we are able to solve these prob- 
lems. 

The Reynolds analogy, first postulated by Osborne Reynolds in 
1890 and verified by experiment, states that momentum and energy 
are transferred in the same way in turbulent shear flow. This anal- 
ogy can he expressed quantitatively by the equation: 



where is the heat flow rate per unit area normal to the surface 
and T 0 is the shear stress at the wall. For fully developed turbu- 
lent flow, the velocity profile normal to the surface is nearly flat; 
hence, integrating between the wall where u is zero and a station 
where u and T can be considered to be mean values 




_Z. 

p du -/dij 



( 2 . 1 ) 



p c r aT /^ 







o Cp(T-T^ a 



r, 2 



{ 2 . 2 } 



or 




( 2 . 3 ) 
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where and are the heat transfer coefficient (Stanton number) 
and friction coefficient respectively. 

If we further write the equations of momentum and energy, 
averaged with respect to time, for mean flcrw parallel to a surface we 



where jl is the viscosity and k i3 the conductivity, and upon re- 
arrangement to better show the similarity: 



number. Thus, it can be seen that if the Reynolds analogy holds for 
the turbulent properties, and if cr = 1, Equation 2. 3 holds through- 
out turbulent and laminar regions. 

For most gases (J~~ — > 0. 7-1.0, so the Reynolds analogy can 
be assumed to hold reasonably well. The analysis and examples be- 
low arc restricted to gas -cooled reactors with turbulent flow in the 



find: 



r 0 = - p uv +/ u <*<y d y 



(2. 4) 




(2. 5) 




( 2 . 6 ) 




coolant channels. 
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2. 2 T emperature Distribution 

Due to non-uniformity of the heat generation in the core of the 
reactor, the average temperature rise of the coolant through the 
reactor core muet be determined to find the overall enthalpy rise 
through the reactor. The temperature rise in the individual channels 
must be found first. Integration of each coolant channel temperature 
rise over the reactor cross-section gives the total enthalpy rise. 

The heat balance for a length of coolant passage, dx = / d£, at 
a given radial station, YJ , is given by: 

rh c p AdT « (2.7) 

where rh is the mass flow rate per unit area, is the specific 

heat at constant pressure, A is the cross-sectional area of the 
coolant passage, T is the static temperature of the coolant, q Q is 
the power density at the point of maximum heat generation in the 
reactor core (usually at the center of the core), £ = is the dimen- 
sionless coolant passage length measured from the channel entrance, 

* —is the dimensionless radius of the core measured from the cen- 
ter, f(£) is the power distribution function in the axial direction, g( ^ ) 
is the function in the radial direction, and ci is the porosity, or void 
fraction, defined as the ratio of the total coolant passage cross-section- 
al area to the total cross-sectional area. In general, c< could be a 
function of Y\ . 

Rewrite Equation 2. 7 as: 

. AXf£> 



(2. 7a) 



I 










-9- 



where 



X- _s 



t 



C f>To 



111 1 

o4 



is a dimensionless parameter. Since m 



is constant for a given coolant passage. Equation 2. 7a can be integrated 
to: 



l 



= T 0 | i + X (?) ) ( f(5) d 5 



( 2 . 8 ) 



The heat balance between the v/all and the coolant is given by 
the relation: 

m c p A iT = C K f u. c p (T w -TJ Pf d $ ( 2 . 9 ) 

where is the Stanton number and depends chiefly upon the Rey- 
nolds number for gases with Prandtl number close to unity; p is the 
coolant density, u is the velocity of the coolant, T is the channel 
surface temperature, and P is the wetted perimeter. 

An empirical relation used in flow through pipes of circular 
cross-section is: 

_ o.oi3 

K —o.t a.U (2. 10) 

cr 

where R e is the Reynolds number. The Reynolds number is de- 
fined as: 

= 4^4 ( 2 . 1 !) 

u~ r 

/ 

where 4A/P = D, , the hydraulic diameter. Since R for a constant 
n e 

area channel varies only with ju. , which in turn depends only on the 
temperature and is approximately proportional to the three- fourths 
power of the absolute temperature, the variation of (and of , 
the friction coefficient) along the passage will be quite small, partic- 
ularly for turbulent flow. 



1 

j 
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Since m a pu , Equation 2. 9 cam be written as: 
dT _ a (t -T) 

' P 1 ~ ; 



( 2 . 12 ) 



where P 3 C h"7T , is constant, and is a dimensionless parameter 
proportional to the heat transfer coefficient, Cj . Hence, the wall 
temperature is given by: 

T w = T + t % (2 - U) 



or, from Equation 2. 8: 

T w = T e +AT - (~40s) * (2 * 14) 

t 

(i T\V 

The maximum wall temperature is given by the condition, — = o, 

1. e. , 

* r 3 ^ = ° ( 2 - 14a > 

From this equation, the value of £ for the maximum wall temperature 

can be found, and substituting into Equation 2. 14 T , the maximum 

m 

wall temperature^ an be obtained. 

2. 3 Pressure Distribution 

The coolant will undergo a pressure drop as it flows through 
the channels. This pressure drop arises from two factors; the ac- 
celeration of the coolant as it passes through the channels, and the 
friction drag opposing the flow of the coolant through the passages. 

The acceleration of the flow results from decreasing density of the 
coolant, which in turn results from increasing temperature and de- 
creasing pressure. The friction drag cannot be avoided since the 
Reynolds analogy holds. 
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The pressure drop along an individual channel can be determined 
by applying the equations of momentum and continuity. The momentum 
balance for an element of length of the coolant passage is given by: 

^ - ? a - 1 c f x =o < 2 - 15 ) 



where p is the static pressure of the coolant and i6 an empiri- 
cal friction coefficient. From the continuity and perfect gas relations: 



<dij_ _ 
' 



met (X) = 



no 



*1* 



(jr) - ?-*■<-'? 



c ij 



-X 4* I 



16) 



where R is the gas constant and subscripts, o, refer to the entrance 
station of the coolant channel. Substituting Equation ^ ^ into Equa- 
tion 2. 15 and rearranging: 




P.M T-| Jg. 
RT.‘ P j 





(2. 17) 



where (3^ = -Cj— , a dimensionless parameter proportional to C^., 

the friction coefficient. 

The second term in the coefficient of dp/d£j in Equation 2. 17 
is a result of the acceleration due to the pressure drop in the channel. 
In many applications this term can be neglected and then Equation 2 ' ^ 
can be integrated exactly; however, in flow through a reactor for pro- 
pulsion applications, the pressure drop may be quite large, so this 
term must be retained. The term containing dT/d£ is a result of 
acceleration due to the heat input to the coolant as it flows through 
the passage, and will be quite large for nuclear reactors. Equation 
X I”? ia an exact equation for the pressure drop, but it cannot be 
integrated in closed form. However, the equation can be integrated 



by an iteration process if the Mach number is low. 



2. 3. 1 Ite ration P rocedure 

Introduce a dimensionless parameter, <£ 

6 = (rn) Z _ (p.u 0 ) RT. _ Q 

Po X Cp.RT.)* RT 0 



defined as: 

« I 



The subscripts, o, refer to the channel entrance, and since the 
product RT is proportional to the square of the local speed of 
sound for perfect gases, €. is seen to be -yM , and thus propor- 
tional to the Mach number at the channel entrance. It will be as- 
& 2 . 

Burned that lL '— = (-Ct\ Jjue << I as well. 

RT, v P> / T. 

Equation 2. 17 becomes, upon substitution of e , 



1" — E- - £ 1 »- "1 fi. , P \ £ _L_ | dT , (3 r 

L P- p/p. 1 Js'-LJ <" £ t. I ag + T ) 

Let p be expanded in a power series in e : 



(2. 18) 



P * + 



f«i JL ^ 

f P - 



{2. 19) 



Substituting this series into Equation 2. 18, we obtain: 

[' + J + J ,2 - 20 > 



+ 4 -'Si*i* T ) = o 



T. 



and , upon collecting powers of €. , the differential equation for p 



(1) 



is found as: 
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and integrating Equation 2. 21, 

"pT = “ Y + 



,(2) 



For the second approximation, the equation for p' ' in: 

A / £ ) -*- A _ i_ ) A. (j!l) 

. p. / 1 p. T # H5VpT; 



( 2 . 22 ) 



(2. 23) 



and integrating Equation 2. 23: 

rx, 



H 

p- 






ft) 



- r , a 



(M 






(2. 24) 



Define the static pressure loss as: 



r fl> * '*> , , \ 

Ap - S - IN + + p ' ' 



C< I 



where pi is the first approximation for p at station 1, the reac- 



<•») 



tor exit, and p v is the second approximation at this station. There- 
fore, the pressure loss ratio may be expressed as: 

&£- = - 6 - «»_£' 

P- P- i TC 

= (2.25) 



or 



AjL. - £• Fj r a. ■*■**“• 

P» 



( 2 . 26 ) 



These functions, Fj and F-, » may be expressed in a more con- 
venient form by integration and rearrangement to the following: 
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(w%)d% + % [(•-?) +(§) d ^_ 



(2. 27) 



TJ 



■J,(r ,-T.) - - [sf ^ 5 



v T(T - T -> J5 






i f 



T. 



-5 . 






(2. 28) 



2. 4 First Order Approximation 

The function, f{£), has been defined as the axial power density 
distribution. Hence, Fj may be written as: 

F . = j? f + A [ B 0 * % ( - B, )] (2. 29) 



where 



B and 
o 



Bj are defined as: 
£f(5)df = const-. 



(2. 30) 



B, -- const*. 

o 



(2.31) 



If we further require that £(§) be symmetrical about the center of 
the reactor, $ » integration of Bj shows that Bj = *B o . There- 
fore, Fj , under these assumptions becomes: 

F, = * X B. ( I + (*/_ ) (2. 32) 



In order to integrate Equation 2. 22 to find the first order ap- 
proximation to the static pressure loss across the reactor. Equation 
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2. 7 can be rewritten as: 

45 - at. -f (§) = k c siv 



<*% 

where K is defined as: 






K 



\T«= 

4,° ^ ^ ~0 

<f p, / R r e 



h$) 



(2. 33) 



and the parameter, X , becomes: 



A -- K 

oc v e 

The exit temperature of the working fluid as a function of radial 
position can be determined by integrating Equation 2. 33, since the 
integration is in the axial direction only: 

44) -T. - KT 0 B, (2- so 

Using these results. Equation 2. 26 for the first approximation to the 
pressure loss can be expressed in either of the following equations: 



^ f t 

P° 






ft - + 

oC ' 



{2. 35) 



or, alternatively, using the relation of Equation 2. 34: 

Ae = fy + ( 1 + fosi) e * -• (2. 36) 

1 0 r 0 

Neither of the two Equations 2. 35 and 2. 36 is in a convenient form 
for integration. However, without changing the value of the equation, 

-j- *y 

we may add and subtract the term, «= — - — C * + , from the 

right side of the equation. — i- 2 - is at present unknown and is 

to be determined. Performing this operation. Equation 2. 35 becomes 

= ^ [P f ^ + 0 + /VJ 6 °4 ^ (2. 37) 
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For a given reactor under given flow conditions, 



T k -T„ 



will be a 



constant. Further, when Integrated over the radial area of the reac- 
tor, Tj and will not differ to a great degree, so the second 
term on the right side of Equation 2. 37 will now be very small when 
compared to either of the other two terms. 

For convenience, rewrite Equation 2. 37 as: 

- E(v^+- [fsfe - D 

P° 



or, rearranging in terms of \/6 : 



- O 



(2. 38) 



(2. 39) 



where: 



|3 f + ( i + ft+A.) 

- (i t/V.lKBj'P’ 



t = 

c 

- 0 +■ (WJj 



Ti - T 

T 



- 



(2. 40) 



Since the middle term in Equation 2. 39 has been purposely made 
very small in comparison with the other two terms, this equation can 
now be solved by successive approximations. Write Equation 2. 39 as 



1 1 rv\ 



(a)\ 






= o 



(2.41) 



and assume a solution of the form; 



- l ,rVV ^ )^| 

I 



(2. 42) 



.After substituting this assumed solution into Equation 2. 41, col- 
lecting similar powers of <5" , and taking the limit as S —>• I » we 
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obtain: 




(2. 43) 



It is more compact to leave Equation 2. 43 in terms of C, D, and E, 
since they are defined in Equation 2. 40. For integration in the radial 
direction, it can be .seen that C contains the variable of integration, 
g( ). For the first approximation, use the first term (which does 
not contain C or D); similarly, for the second approximation to fe. 
terms linear in C and D only are required. Higher order requires 
higher powers. 

2. 5 Average Temperature Rise Through Rea c tor 

The total enthalpy rate of increase, AQ , in the exhausting gases 
is: 



where rh is defined as the total mass flow rate found by integrating 
the individual coolant mass flow rates over the cross-section of the 




(2. 44) 




However, AQ may be found directly as: 
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Aq = 



TV* <- P j ' '«j ^ d *| 



i-TTo 2, ( 8, J ^ 



(2. 45) 



Therefore, by equating Equations 2. 44 and 2. 45 we may 3olve for the 



average temperature rise through the reactor: 



T t ~T 0 = 






•p [ ©< (>] ) > }^ > ] 



(2. 46) 



- KT.B .-4^ 

J^f^jy/e 7>j) 

The exit temperature at the centerline of the reactor io found 
by putting rj equal to zero in Equation 2, 34: 



TM-r. - Kr.JU^. 



oc 



b>. 



(2.47) 



- A (o'! r 0 b. 



(2. 474) 



Hence, dividing Equation 2. 47 by 2. 46, the ratio between the center- 
line exit temperature of the reactor and the average exit temperature 

ia: /— ■ 

/ g n£y | 



( o) - T 0 
T, - T„ 



$<«> 



<* /efoj' j rj 



(2. 48) 



It is not convenient to express the relation between the maximum 



wall temperature and the average temperature increase across the 
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reactor in general. However, if f(£) and g{ rj ) 
this relation can be found quite readily as will be 
examples. 



are known functions, 
shown in the following 



; 
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in. EXAMPLES 



3. 1 Reactor Calculation Using Spe c ific Power Dist rib ution Functions 

Thus far in this thesis the general equations for determining the 
overall temperature rise and pressure drop through the reactor have 
been obtained while leaving the specific power distribution as general 
functions. An example will now be given showing how the results ob- 
tained in the general derivation might be applied to a given reactor. 

The general solutions for the neutron density distribution in a 
cylindrical nuclear reactor are given in terms of trigonometric and 
Bessel functions. However, for convenience, a parabolic neutron 
density was assumed since the integrated area under a parabola dif- 
fers but slightly from the more generally accepted solutions of Bessel 
functions or trigonometric functions. For this example, aesume a 
right- circular cylindrical reactor having parabolic normalized power 
distribution characteristics as: 



Further, assume the porosity is uniform across the reactor. The 
power distribution parameters to be perturbated will be & and 
which determine the relative "flatness" of the axial and radial pov/er 
input cf the reactor. These parameters are primarily a function of 
the amount of refloctor material used in the reactor design if homo- 
geneous fuel loading is assumed. 




(3. 1) 



<?<>?> * 1 



(3. 2) 
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With these functions, B is found from Equation 2. 30 to be: 



B„ = 5 = I - % 



(3.3) 



Using Equation 2. 48 and working out the integrals indicated to 
the second approximation in /dP and with: 

% (o) ~~ 1 



y/6(0) = /A P /p 0 



( I 









(3.4) 



we find, after simplification, that: 



T» -T. 
T~T a 



f 



i - \ 



t,-t. , 

ft * 



ft 



r(3. 5) 



To find the maximum wall temperature, which will occur on the 

centerline of the reactor, use Equation 2. 14a to determine the value 

of £ to put in Equation 2. 14. From Equation 2. 14a, T occurs 

m 



when 



s-i-t-jf 

and hence, from Equation 2. 14: 



l 

46 



r ”s±— . XU) 

L 



/i - 



1 



i _ (9 
Z 



(3.6) 



and using Equation 2. 47a 



- l« 



T, (o) - i 0 1 e 4 



^ + ^ ^5(3 



i e 

'•a 



(3.7) 
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Hence, given a maximum wall temperature limitation in the 
reactor, the relation (as a function of 0, <(>, p, andp^) between the 
average temperature increase of the gas flowing through the reactor 
may be found by multiplying Equations 3. 7 and 3. 5: 



- ip, 



r 



0 - %i 0 - &'■ 



1 



|9 f » 








(3.8) 



It is simpler to leave Equation 3. 8 in this form; however, for given 

values of the parameters of the e ^nation, it reduces to a simple 

T — T t _ T" 

cuadratic to determine — 1 — — for a given — — — ■ 

G X 

If all wetted surfaces of the channel are heat transfer Ftrfaces, 
and from the Reynolds analogy, it is reasonable to assume that p » 0^ 
for fully developed turbulent flow. Furthermore, it is assumed that 
A = 8 for comparison purposes, though neither of these assumptions 
are requirements of the analysis. The results obtained from using 
Equation 3. 8 are shown in Table I. This table show s the maximum 
wall temperature required to obtain a temperature increase of 
1000°R and 1500°R for an initial gas temperature of 1000°R for var- 
ious values of p and <{>. The tabic was prepared on the assumption 
that <j> s 8 and p = p^ . The required temperature increases were 
selected quite arbitrarily, but they furnish an indication of the trend 
in maximum wall temperature versus average temperature increase 
through the reactor. From Table I, it can be seen that the maximum 
wall temperature required for a given temperature increase of the 
working gas is markedly increased as the power distribution varies 
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av/ay from a linear relation ( <*> equals zero) or when the coefficient 
of heat, p , gets very low. 

3. 1. 1 No -orifice Condition 

The static pressure loss may be found by using Equations 2. 44 
and Z. 45: 

YKc p (VT 0 )= aq 

= 2 rra* f>„ /rt; KT 0 b - JJ t 1 ] ) dnr \ 

(3.9) 

^ -T: ;.~J» ^ p„ \/rt; K (> - **/>) Ci - */*.) 

t ra. x T 0 



But, from the simplification procedure (not shown) used in deter- 
mining Equation 3. 5, it was found that 



K 



pC T i~ T ° / __ A P/P° 

0 - %) ( ' - *4) v 1- -^H 1 */^/0 

* o 



Therefore, Equation 3. 9 can be expressed as: 



AE _ J $ + _h Tz— ( i •<- u/\Lf yy^ \ 



T". 









(3. 10) 



Yx 

However, — is just the mass flow rate, p u , per channel, 

o. TT 0 J~ * r O O r 

bo Equation 3. 10 can be reduced further to: 

A fL ^ J 1, , ( 



■ 3 t ’ 



I +• 



r. 



*/-n e 



(3. 10a) 



\T& 



Figure 1 shows the variation in static pressure loss with 

(which i3 proportional to the mass flow) for various values of 



for an assumed temperature rise ratio of unity across the reactor. 
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3. 1. 2 Artificial Flow Restriction (Orifices ) 

For installed orifices scaled so that the macs flow in the outer 
channels is restricted to the correct value to give a uniform tempera- 
ture rise across the reactor, the average exit temperature will be 
identical with the exit temperature from the center channel of the un- 
restricted reactor. Therefore, from Equation 3. 7 

1 o l 



T. 



-T. "-“Ht 



(%V iWi 



H'- la 






2. (a 



(3. 11) 



If the orifices restrict the flow correctly, the pressure decre- 
ment they provide will make all channel pressure losses identical 
with the pressure loss In the centerline channel, which is unrestrict- 
ed. Therefore, the first approximation to the pressure loss ratio 
can be determined from Equation 3. 4 by assigning a value to \fe(d) 
under given reactor conditions. 

The other change that will result from artificially restricting 
the channel entrances will be the total mass flow. This total mass 
flow can be determined directly using Equation 3. 11. Since the 
average temperature of the gases is the exit temperature, a constant 
across the reactor, by using Equation 2. 34 we are able to determine 
m if it is recalled that it is a factor in the parameter, \f&' . Hence, 
integrating, the total mass flow is : 

rq - zfr ^J J ^ >)*<■) 



T* 



yq = 4 



l 



C pT* 



T, 



T„ 






(3. 12) 



However, since T, 



r« s 



is tlie exit temperature of the centerline chan 



nel, Equation 3. 12 can be rewritten by using Equation 2. 47 and the 
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fact that m 

Yn 




- oiTTcf-^ 



\f £(°) ( 1 ~ ^2.) 



(3. 13) 



3. 2 Rocket Example 

Using the assumed reactor power distribution characteristics 
of Equations 3. 1 and 3. 2, a rocket engine using molecular hydrogen 
as the working fluid was analyzed. No dissociation effects were con- 
sidered. The fuel was assumed to be liquified in the propellant tanks; 
however, due to cooling requirements of the nozzle and reactor reflec- 
tor, the hydrogen was assumed to have reached a gaseous state before 
entering the pumps of the "combustion chamber". 

Account was taken of the pov/er required to drive the pumps by 
using a portion of the gas from the reactor exit to drive a turbine 
which furnished power for the pumps. The temperature of the gases 
as they were bled from the reactor exit was reduced to acceptable 
turbine operating temperatures by mixing it with the required amount 
of coolant hydrogen to lower the turbine gas inlet temperature to an 




,o; <0 (z) 

Schematic Drawing of /foe kef 
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For this example, define the following parameters: 
nHp mass flov/ rate of pump 
mass flow rate of nozzle 

<5" bleed fraction to drive auxiliary turbine, therefore, 

nhp = ( i + <T) 

7T^ reactor exit-entry total pressure ratio = p iT /p IT 

Tr m assumed exit-entry total pressure ratio due to mbdng 
loss from exit temperature gradient = pz T /p x 

"TC"r\ nozzle exit-entry total pressure ratio (assumed 1. 0) 



The thrust, F, of a rocket is given by the equation: 

F = rrV ^4 + ^A ( Pa ~ Pa-l 

+ 4 r (3. 14) 

where p a is the ambient pressure and subscript 4 refers to e?dt 
station. The specific impulse, I , is found by the equation: 

I - - -ili- , ( P4 - M 

9 f "V At . < 3 - 15 > 



Using the defined parameters, the nozzle mass flow rate, , 
can be found using the continuity equation and the perfect gas relations 
for a choked throat as: 



rw 



■f _ 










$ -H 
l-(ff-i') 



(3. 16) 



If the expansion ratio, 4> 4 t is specified, the nozzle -exit 
ratio can be determined from the following equation if isentropic 
flov/ is assumed in the nozzle: 



(3. 17) 




where ^* r /p A ~ 11 00 ^ Xt /P+ • Equation 3. 17 uniquely defines the 



area ratio and subsequently the exit Mach number from the relation: 

S + » 

A, 






6 1 1 



* - KA 

z 






*(«-•) 



(3. 18) 



Kence, is now known. 

Since has now been determined, the exit velocity, , can 



be found if is known. 



T < - T v ( 



, + 

z 






But, 



T»t - t 



Therefore; 



t 4i - t, t - T, t ( > 

-- t, t fi * 

L [ { T 

t 4 T * t .(' - *?)( • - 



(3. 19) 



where lj ~ - — is the static temperature increase found from Equa- 

' l 

ticn 3. 8, T , is assumed equal to T . • since the entering Mach num- 

i I rjt 

ber will be very small, and was determined from Fig. 6, Dailey 
and Wood (2), modified to account for the friction in the coolant pas- 
sages. 

Knowing is now known, and the exit velocity, u^, is: 



U ~ U 4 M 4 * U 



(3. 20) 
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The total pressure ratio across the reactor is found in a similar 



manner using the perfect gas relations 

Ap r 



7T, 



R 



I 



P, 



3 + 



i « -i 



{3. 21) 



where 



Ap 



is the static pressure ratio across the reactor as deter- 



mined from Equation 3. 10a by assigning various values of <£ . After 
having determined the static pressure of the exhaust is given by 



the relation: 



ft = P4 T (‘ * ■ ± i Lt K) 



& 

*-( 






J - 1 



(3. 22) 



The total mass flow rate that must be pumped by the pumps is 
known if e. is assigned, since from Equations 3. 10 and 3. 10a: 

ty i p 



TTa- 



P' 



£ rt; 



(3. 23) 



In order to pump the required total amount of mass flow, the pumps 
expend power which was assumed to be supplied by the turbine. Equa- 
ting the turbine power output to the required pumping power: 



Pit 

p > 



<Zy ^ ^turb . 



(3. 24) 



Therefore, assuming that the mixture of bled gas and coolant gas 
enters the turbine at 2000°R and exits at the reactor entrance tem- 
perature, the amount of bled gas from the reactor exit is given by: 



rr\ n - 



P'T Ht 



(3. 25) 
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in working out the example, a "rubber” engine was assumed, 
i. c. , the rocket was assumed to be correctly designed for each vari- 
ation in the parameters. A. graphitic reactor is used as the basis for 
assuming the high maximum wall temperature (Ref. 3) since it was 
visualized that a rocket engine would not have a long operating life. 

No account was taken of any possible mass addition to the flow from 
the reactor material itself. 

Assumed conditions 

Pj = 50 atmospheres 

p = 0. 1 atmospheres ( 53, 000 feet) 

3i 

T = 5500°R 

w 

m 

tt = 0. 90 for no-orifice condition 

IT. 

1. 00 for orifice condition 
a 0. 5 

Tj = 500°R 

F V*- so 

fra 1. 35 

C = 3.3 Btu/lb °R 

P 

Using these conditions, the exit-throat area ratio is 5. 53 and the 
exit Mach number is 3. 17. 

Figure 2 shows the specific impulse versus p obtained from the 
calculations for 4> = 0 = 0. 5. A linear power distribution (<J> = zero) 
is also shown for comparison. No attempt is made to distinguish 
between the values of mass flow used at each p , since their effect on 
the specific impulses was negligible, although their values do affect 
the pressure loss across the reactor. From these results, the spe- 
cific impulse appears to be primarily a function of the temperature 



-3‘ > - 



rise imparted to the gas in passing through the reactor. Therefore, 
to get high specific impulses in a rocket, it would seem preferable to 
attempt to increase the temperature of the coolant gas by increasing 
the heat transfer coefficient, although this lowers the allowable mass 
flow rate to preclude very high pressure losses, or choked flow, in 
the coolant passages due to the higher heat input to the gases. .For 
this same reason, orificing appears to be advantageous in increasing 
the specific impulses. 

3. 3 Ramjet Example 

In the ramjet example, the conventional burner is replaced by a 
nuclear reactor having the power distribution characteristics of Equa- 
tions 3. 1 and 3. A This analysis was based on standard ICAO condi- 
tions at 1000 feet altitude. 

In working out the example, a "rubber" engine was assumed and 
the nozzle was assumed to be designed to expand to atmospheric pres- 
sure unless the throat static pressure was lower than ambient. In 
such a case, no expansion was considered, and the nozzle was assumed 
to be terminated at the throat. 
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(Z>) (s') 



SchcnnaTic. f)raw/ny o-f 



(t) 



(s') 



-31- 



For tliis analysis, define the following parameters: 






tr r 



ratio of free -stream total pressure to static = F^r/p 0 
ratio of total pressure across diffuser » f’a.-r/p 
ratio of total pressure across Reactor = f *r/p tr 



ratio of total pressures due to assumed mixing loss 

(sumed 1 for 



in non- restricted case = PsV/p (Ass 
restricted case. ) ' 3r 

ratio of nozzle total pressures (assumed 1) 



The thrust, F, of a ramjet is given by: 

f = - «. -f) 1- | - p./ ps ) 



(3. 26) 



Due to the fact that ramjet engines must operate for long peri- 
ods of time, and in an oxidizing atmosphere, a different type reactor 
material must be used in a ramjet than in a rocket. Because of this, 
the materials suitable for ramjet applications can withstand a much 
lower maximum wall temperature than those for a rocket engine. The 
heat input to the coolant gases will be much less, and consequently the 
acceleration of the gases in the coolant passages will be less. There- 
fore, in this example, it was assumed that the static temperature 
rise and the static pressure loss across the reactor were approximately 
equal to their total, respectively. 

can be found by finding the ratios of exit Mach number to 
the free stream Mach number and exit total temperature to the initial 
total temperature. 




(3. 27) 



2 - 



But: 



JSu 

C't 



f>s r -Mi * h ** "■* ^ (' + 4 M ") < 3 - 28 > 



and from this we obtain: 



K 



^ It 
. M. I 



- 

Ps 



i \ + Sr i 



H 1 M « l 



t! M *■ 



(3. 29) 



Also: 



T c 



•T _ 



T, 



= \ + 



AT 

T_ 



(3. 30) 



where ^=f = — 1 _J° - of the general reactor analysis. 

°t *• 

Therefore, the exit velocity can be expressed as: 



u< 



F 



l + 



AT 

To, 



r 



Ld M l 
2 - ° 






Since it was assumed that the static pressure loss and total 
pressure lo3S are approximately equal, ir^ is found by using Equa- 
tion 3. 10a to determine Ap/p for various values of £ . 



^ - 1 - At> 4 M 



(3. 32) 



Knowing - , the throat pressure can be found from: 



h = 

'?» v fii-t/ 



K K, ^ 



(3. 33) 



If this ratio is less than unity, no further expansion was assumed; 
if the ratio was greater tlian unity, ideal expansion to ambient pressure 
was as sximed. 
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Figures 3 and 4 3 hov/ the results obtained for C_ (based on the 

i? 

reactor cross-sectional area), tho thrust coefficient, versus Mach 
number for various values of p at an assumed maximum reactor wall 
temperature of 2960°R. Cp is defined as the total thrust divided by 
the free stream dynamic pressure and the cross -sectional area of the 
reactor. Values of used in the example were: 



M 0 




2 S 


o.8s 


2.S — S' 


?£ 


3.S — 4. S' 


a ® £ 



tr_ , the assumed mixing loss, was arbitrarily assigned a value of 
0. 90 for the non- restricted example and unity for the orificod exam- 
ple. Results shown are for an assumed reactor porosity of 0. 5. 

The curves shown are for £ = 0. 03, but other moss flow 
rates give similar curves. The one shown is about optimum, since 
it was found that further increase in mass flow rate increased the 
pressure loss across the reactor to a point where tho overall thrust 
is reduced. While it should be stressed that the initial conditions and 
reactor configuration chosen were arbitrary and this example was 
worked to demonstrate the methods, still, the trend of the results 
should be realistic, if not their absolute magnitudes. From these 
two figures at the arbitrarily chosen configuration, optimum per- 
formance is obtained at j3 of approximately 3-5 at a Mach number of 
about 3. Figure 4 shows that v/hile Cp at the optimum Mach number 
is slightly increased by orificing, the vehicle performance drops 
more sharply away from the optimum Mach number, so it is doubtful 
that there is much advantage to orificing for a ramjet. 
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TABLE I 



Maximum Wall Temperature Required to Obtain Given Average 

Temperature Increase of 

Coolant Gas (Initial gas temperature assumed to be 1000°R) 





AT 


= 1000 




AT = 1500 














P -= P £ 


Tj(o) 


T 

w 

m 


Tj(o) 


T 

V/ 

m 




20 


3430 


3440 


4870 


4880 




10 


3460 


3470 


4930 


4940 


4 » = l. o 


5 


3520 


3625 


5040 


5250 




1 


3860 


7000 


5600 


10, 630 




20 


2790 


2818 


3770 


3320 




10 


2801 


2868 


3810 


3910 


<*> = . 75 


5 


28 28 


2981 


3840 


4080 




1 


2950 


4760 


4030 


6850 




20 


2418 


2518 


3160 


3310 




10 


2420 


2533 


3170 


3340 


$ a . 50 


5 


2430 


2620 


3190 


3440 




1 


2480 


3720 


3260 


5170 




20 


2173 


2450 


2771 


3190 




10 


2175 


2452 


2772 


3190 


<J> = . 25 


5 


2177 


2490 


2779 


3250 




1 


2193 


3190 


2801 


4315 
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TABLE I 
(continued) 



AT = 1000 


AT = 1500 


0 = P f 


Tj(o) 


T 

w 

m 


Tj(o) 


T 

V/ 

m 


20 


2000 


i 

2000 


2500 


2500 


10 


2000 


2000 


2500 


2500 


5 


2000 


2000 


2500 


2500 


1 


2000 


2000 


2500 


2500 



airo 2 
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FIG. I - STATIC PRESSURE LOSS VS. DIMENSIONLESS 

MASS FLOW RATE FOR UN R ESTRICTED COOLANT 

PASSAGES 
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FIG. 2 - SPECIFIC IMPULSE VS. HEAT COEFFICIENT FOR 
NUCLEAR - HYDROGEN ROCKET EXAMPLE 
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0 2.0 3.0 Mq 4.0 5.0 



FIG. 3 - RAMJET EXAMPLE (UNRESTRICTED COOLANT 
PASSAGES) - THRUST COEFFICIENT VS. FREE 
STREAM MACH NUMBER 




FIG 4- RAMJET EXAMPLE (ORIFICED COOLANT 
PASSAGES )- THRUST COEFFICIENT VS. 
FREE STREAM MACH NUMBER 
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